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We numerically investigate the conical emission (CE) from ultrashort laser 
filaments, both considering and disregarding the higher-order Kerr effect 
(HOKE). While the consideration of HOKE has almost no influence on the 
predicted CE from collimated beams, diflFerences arise for tightly focused 
beams. This difference is attributed to the different relative contributions 
of the non-linear focus and of the modulational instability over the whole 
filament length, (c) 2013 Optical Society of America 
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The filamentation of ultrashort laser pulses attracts much interest due to its spectacular 
potential applications [1] such as remote sensing [2], lightning control [3], water condensation 
[4], or the generation of THz radiation [5]. A prominent property of filaments is conical 
emission (CE) [6,7], which denotes the spatio-spectral distribution of the supercontinuum 
generated in the filaments, with a well-defined relationship between the wavelength and the 
emission angle. It provides an angle-resolved spectrum refiecting the complex dynamics of 
femtosecond pulses propagating in nonlinear dispersive media [8]. CE is observed on both 
sides of the fundamental wavelength [9]. This spectral symmetry is understood by considering 
either modulational instability [10,11] or the X-wave propagation regime [8]. 

Filamentation is generally described as a dynamic balance between Kerr self-focusing 
and defocusing by the plasma generated at the non-linear focus [1, 12-14]. The temporally 
asymmetric distribution of the plasma along the pulse local time is expected to generate a 
spectral asymmetry of CE in the region of the filament onset (or beam collapse). [7, 13-15]. 

However, a controversy is currently running about the possible contribution of higher- 
order Kerr effect (HOKE) [16] to filamentation in gases [17-20]. Since their consideration 
leads to expect 10 to 100 times lower plasma densities in filaments [17], and CE is both a 
key property of filaments and expectedly strongly connected to ionization, it is important to 
understand how the HOKE would impact CE. 
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In this Letter, we numerically investigate both the impact of beam focusing on CE and 
the build-up of the CE along the filament length. We find that both models reproduce CE 
on both sides of the fundamental wavelength for collimated beams, but that considering the 
HOKE is necessary to reproduce the observed CE from focused beams on the high-frequency 
side of the spectrum [21]. 

We simulated the propagation of a 5 mJ, 45 fs pulse (13 critical powers Per) centered at 
800 nm, launched collimated in air with a radius of 3.8 mm at 1/e^. The code is based on the 
unidirectional pulse propagation equation (UPPE) [22] , in which we used the published values 
for the HOKE [16]. The model, described in detail in [18] in the case of argon, was completed 
by a detailed weak-field quantum modeling of the Raman-induced rotational eflFects in the 
air [23]. The far-field CE pattern is computed as Fourier (in the temporal domain) and 
Hankel (in the spatial domain) transforms of the electric field. 

Both the behavior of the conical emission and the impact of the HOKE consideration 
strongly depend on external focusing. For a collimated beam, both models reproduce the CE 
from a collimated beam at the experimentally measured angles on both the visible [24] and 
the infrared [9] sides of the spectrum. However, predictions diverge when the focusing gets 
tighter. The truncated model predicts the disappearance of the CE, first on the visible side 
(Figure Id), then on the infrared side (Figure If) of the fundamental wavelength. Conversely, 
the full model predicts CE on both sides of the spectrum, regardless of the numerical aper- 
ture. Focusing only offsets the CE angle by an amount comparable to the numerical aperture 
imposed to the beam, and broadens the angular distribution for each spectral component 
(Figure la,c,e). Let us note that CE was indeed observed in the visible side of the spectrum 
in beams with comparable /-numbers by Faccio et al. [15]. Such behavior is only repro- 
duced by the full model, although experimental conditions differ slightly from those of our 
calculations, preventing a direct comparison. 

The similarity of the conical emission patterns predicted by both models from a colli- 
mated beam (Figure 2) apparently contradicts the recent finding that considering the HOKE 
in models suppresses the CE from the filament onset region [20]. This paradox can be re- 
solved by considering the build-up of CE over the propagation distance (Media 1 of Figure 
2 and Figure 3). Note that, in the truncated model outcome, the intensity spikes are not 
numerical artifacts. They stem from refocusing of the the trailing part of the pulse during 
the propagation [25]. Although both models predict the same filamentation onset distance 
(z 9 m. Figure 3a), the truncated one yields a slightly earlier rise of the CE, which is more 
pronounced on the visible side of the spectrum (Figure 3b, c). This is especially the case at 
wavelengths close to the fundamental one, where the build-up of the CE appears offset by 
up to 1 m between the two models (See inset of Figure 3c). The offset in the CE build-up 
predicted by the two models at the filament onset may in principle be used to distinguish 
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between their validity [20]. However, this offset has the same order of magnitude than the 
mismatch displayed the filament onset location observed experimentally and predicted by the 
full model on one side, and its location predicted by the truncated model on the other side 
(See Figure la,b of [20]). We therfore expect that such test is hardly conclusive in practice. 

The differences observed between the two models including and disregarding the HOKE 
can be understood by considering that two processes contribute to the CE, namely the 
modulational instability (MI) all along the filament length, and refraction of the continuum 
by the refractive index gradient induced by the transverse plasma density profile. The gain 
of the MI contribution to CE writes [11]: 



where D = ^ is the spatio-spectral dispersion operator, uj is the frequency detuning 

relative to the fundamental frequency, fco (taken at cjq) and the total and transverse wave 
vectors, k^^^ and fc^^^ the first- and second-order dispersion term, /q the incident intensity 
and n2,efr the effective non-linear refractive index. For n2,eff < 0, the gain cancels and the 
supercontinuum emission is confined on the beam axis [11]. Conversely, since the electron 
density accumulates over the pulse duration, its contribution intrinsically bears a strong 
temporal asymmetry, which converts into a spectral one, because dispersion and self-phase 
modulation red-shift the front of the pulse and blue-shift its trail. Furthermore, the plasma 
contribution can be considered as negligible in the framework of the full model, where the 
plasma density is insufficient to significantly influence the pulse propagation [18]. 

Compared to the full model, the truncated one predicts a nonlinear refractive index in- 
version at a higher intensity [17], resulting in a higher MI gain, hence in an earlier growth 
of CE at the filament onset, followed by a larger intensity. Furthermore, the high-intensity 
peak predicted at z 9 m by the truncated model (Figure 3a) will result in high plasma 
density at the filament onset, offering a supplementary boost to the early CE emission on 
the visible side of the spectrum. 

Tighter focusing results in higher plasma densities, especially in the framework of the 
truncated model. This plasma density accumulating during the pulse results in an increas- 
ingly negative contribution to n2,eff, ultimately inverting it to negative values preventing CE. 
Focusing tighter will yield to an earlier inversion of n2,eff , extinguish CE from MI first on the 
visible side of the fundamental, and then on the infrared side, as displayed in Figure ld,f. 

In conclusion, we have investigated the influence of beam focusing on conical emission as 
well as the dynamics of its build-up along the filamentation of ultrashort laser pulses. CE 
is emitted from the whole filament length through MI [11]. For collimated beams, which 
yield long filaments, the contribution of the non-linear focus region is moderate and the 
consideration of the HOKE has no practical impact on the predicted CE. For tightly focused 
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beams, the high plasma density predicted if the HOKE are disregarded would lead to expect 
the extinction of CE. Our work therefore illustrates the need to consider CE emission from 
the whole filament length. It also suggests the behavior of CE at tight focusing as a successful 
experimental test of the need to consider the HOKE to adequately describe CE from laser 
filaments. 
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and Maioli et al. [24] on the IR and visible sides of the spectrum, respectively. 
Media 1 displays the build-up of these patterns over 8-14 m propagation range. 10 
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Fig. 1. (Color online) Far- field angularly-resolved spectrum of 1 mJ, 30 fs pulses (4 Per) 
centered at 800 nm, with a radius of 3 mm at (a,b) / = 4 m, f/# = 1333; (c,d) / = 2 

m, f/# = 667; (e,f) / = 1 m, f/# = 333, in the case of the full model considering the HOKE 
(a,c,e) and of the truncated model disregarding the HOKE (b,d,f) 
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Fig. 2. (Color online) Far- field angularly-resolved spectrum of 5 mJ, 45 fs pulses centered at 
800 nm, after 24 m of propagation in air, including 16 m of filament at ion. (a) Full model 
considering the HOKE; (b) Truncated model disregarding the HOKE. White dots display 
the experimental data of Theberge et al. [9] and Maioli et al. [24] on the IR and visible sides 
of the spectrum, respectively. Media 1 displays the build-up of these patterns over 8-14 m 
propagation range. 




Fig. 3. (Color online) Build-up of the conical emission from a 5 mJ, 45 fs pulse centered at 
800 nm. (a) Total on-axis intensity, (b-d) Intensity of the CE lobe at (b) 500 nm; (c) 700 
nm and (d) 1500 nm as a function of propagation distance, as predicted by the full model 
considering the HOKE and the truncated model disregarding the HOKE. Note that panels 
b-d share a common intensity scale. 
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